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Abstract 

Tlie accumulation of electrons and positrons in the vacuum magnetosphere of a neu- 
tron star with a surface magnetic field of S ~ 10^^ G is considered. It is shown that 
particles created in the magnetosphere or falling into the magnetosphere from outside un- 
dergo ultra-relativistic oscillations with a frequency of 10 — 100 MHz. These oscillations 
decay due to energy losses to curvature radiation and bremsstrahlung, with their frequen- 
cies reaching 1 — 10 GHz. Simultaneously, the particles undergo regular motion along the 
force-free surface along closed trajectories. This leads to the gradual accumulation of par- 
ticles at the force-free surface and the formation of a fully charge-separated plasma layer 
with a density of the order of the Goldreich- Julian density. The presence of a constant 
source of electron-positron pairs in the magnetosphere due to the absorption of energetic 
cosmic gamma-rays leads to the growth of this layer, bringing about a rapid filling of the 
pulsar magnetosphere with electron-positron plasma if the pair-creation multiplication 
coefficient is sufficiently high. 
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1. INTRODUCTION 

The generation of electron-positron plasma in the magnetosphere of a rotating,magnetized 
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neutron star has been well studied [lH3[ • This has been used as a basis for developing the 
physics of a stationary pulsar magnetosphere filled with plasma whose density appreciably 
exceeds the so-called Goldreich- Julian density [sl, ncj = |f2 ■ B|/27rce, corresponding to 
stationary rotation of the magnetosphere right out to the light cylinder, Rl = c/Q. Here, 
Q is the angular velocity of the star's rotation, c the speed of light, and e the positron 
charge. However, it is difficult to understand what mechanisms are generating the radio 
emission of pulsar magnetospheres, where the plasma is created, based on observations of 
radio pulsars operating in a steady-state regime. Tracing the dynamics of the development 
of the emission at different frequencies could prove important for our understanding of 
physical processes occurring in radio pulsar magnetospheres Q]. 

Moreover, many observations of non- stationary radio pulsars have recently appeared. 
These concern first and foremost so-called switching radio pulsars, whose radio emission 
is observed only during certain time intervals appreciably exceeding the rotational period 
of the star, such as the pulsars PSR B1931+24 ^ and PSR J1832+0029 [6]. In addition 
to switching pulsars, a group of so-called nulling pulsars has long been known, which 
likewise display no radio emission during certain intervals, but with these being not as 
regular as those for switching pulsars 0]. Another group of non-stationary radio sources 
has recently been observed: rotating radio transients (RRATs), which are sporadically 
flaring radio sources. The phases are preserved during these flares, and the corresponding 
measured periods are characteristic of ordinary radio pulsars |8|. There is no doubt 
that these are also rotating neutron stars. In our view, all these sources exhibit a non- 
stationary generation of plasma in neutron-star magnetospheres. It is therefore important 
to understand how the magnetosphere of a rotating, magnetized neutron star is fllled with 
plasma, which is the reason for the operation of radio pulsars. 

Together with jol, the current study is concerned with an initial investigation into 
the "ignition" of the magnetospheres of neutron stars — the dynamics of the flUing of the 
vacuum magnetosphere with electrons and positrons created in the magnetosphere. In ol, 
we analyzed the dynamics of the motion of particles in the vacuum magnetosphere of a 
rotating neutron star possessing a strong magnetic fleld 10^^ G). Here, we consider the 
deceleration, capture, regular motion, and accumulation of particles near the force-free 
surface in the magnetosphere — the surface where the electric fleld along the magnetic-fleld 
lines vanishes. 



2. DECELERATION OF THE PARTICLE MOTION 

In our study of the motion of charged particles j9|] , we showed that particles approaching 
the force-free surface begin to undergo oscillations. During these oscillations near the 
force-free surface, the charged particles lose energy to radiation. Let us consider the 
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time dependence of the rate of decrease of a particle's energy. The particle's energy 
is characterized by the Lorentz factor possessed when the particle crosses the force-free 
surface, which we denote C. Our main interest is the variation of the energy in the 
ultra-relativistic case, when C ^ 1; in the non-relativistic case, C ^ 1, and this quantity 
remains essentially constant. 

We obtain from the Dirac-Lorentz equation 



dC 2 



where 7 is the particle's Lorentz factor, p the radius of curvature of the magnetic-field 
line along which the particle moves, and a = /he ~ 1/137 the fine-structure constant. 
Here and below, we use dimensionless units, measuring the strengths of the electric and 
magnetic fields in units of the so-called critical field i?cr = mlc^/eh ^ 4.4 x 10^'^ G, the 
particle velocity in units of the speed of light c, the particle charge in units of the positron 
charge e, the particle mass in units of the electron mass mg, the particle energy in units 
of the electron rest-mass energy rrieC^, all distances in units of the Compton wavelength 
of the electron X = h/rrieC ~ 3.86 x 10^^^ cm, and all times in units of X/c. Note that, 
in these units, 1000 km ^ 2.6 x 10^^ and 1 s ^ 7.8 x 10^°. Then, C = 7 + uj^{x[)^/2 = 
1 -|- where u is the frequency of the non-relativistic oscillations, x[ the current 

longitudinal coordinate of the particle along the magnetic-field line, and A the amplitude 
of the oscillations. The frequency of the non-relativistic oscillations is given by 

2_ bo ■ V(Eo ■ Bo) 
^ 5 ' v^J 

where Eq, Bo, and bo = B0/-B0 are the electric and magnetic vectors and a unit vector 
along the magnetic-field line at the point where it intersects the force-free surface. To 
order of magnitude, the frequency of the non-relativistic oscillations is w ~ ^^/B/Rl, 
which corresponds to the frequency u = 00/2?: ~ 1—10 GHz. 

If we neglect the radiation-reaction force, which formally corresponds to a = 0, then 
C = const. In this case, C will be an exact first integral. However, generally speaking, 
there is a radiative reaction, and the particle will lose energy. This leads to a decrease in 
C. The rate of decrease of C will not be constant over times of the order of the period 
of the particle motion, because it depends on the specific Lorentz factor of the particle 7, 
whose value must obviously change over one-quarter of the oscillation period from 1 at 
the turning point where x[ = A to C ^ 1 at the point x[ = 0. However, this will not 
hinder our investigation of the variation of the particle energy, as long as C itself varies 
little over the oscillation period: 

dC C 

Satisfaction of this condition of applicability of the adiabatic approximation means that we 
can accurately average ([1]) over the oscillation period. We must average expressions that 
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depend on the particle coordinate x\ not formally, but instead using the explicit relation 
x\{i) = A{—1)^^^^ (z — 2h(z)) for the time dependence of the coordinate of the oscillating 
particle, which is exact when C ^ 1. The function h{z) is given by h{z) = \_{z + 1)/2J, 
where z = t/A and [y\ denotes the integer part of y. 

Let us suppose we have a function /, which depends on the distance from the equi- 
librium point, with the remaining parameters being arbitrary, with the only specification 
being / = f{\x[{t)\). It can then easily be shown that an average over the period T can 
be replaced by an average over one-quarter of the period A. If we express / as a function 
of z = t/A using the transform / = f{z) = f {\x[{Az)\) , the mean of / over the oscillation 
period will be given by 

1 

(/) = I f{z) dz. (4) 



The angular brackets denote the mean of a function over the oscillation period. To average 
the right-hand part of ([I]), we must know the Lorentz factor 7: 

^C[l-{z-2h{z)f]. (5) 

Formula ([5]) was obtained using the relation C ~ uj'^A^/2. Further, we neglected unity 
compared to C in the expression for 7, which does not affect the averaging results, since 
C ^ 1. After averaging ([T]) using (jl]), we obtain the differential equation for C 



dC 4 
— — = — auj 6 
dt 9 
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The exact solution of (jH]) has the form 

e{t)=[{l + e,')e'''^^-lY"\ (7) 

where e{t) = C{t)/ Ccurv is the normalized first integral and Eq = e{0) its value at the 
initial time t = 0. Here, we have also introduced 

Ccurv = ^^{Oopr^' ^ 1.2 {Oopr/' (8) 

and the decay time constant 

4 acj^ 

The physical Td are not difficult to understand based on the 

behavior of the solution ([7]) for various initial values Eq. Let first Eq <^1. The solution of 
^ is then given by 

E{t) = E^e-"^'' {Co Ccurv). (10) 

We can see that the case Cq ^ Ccurv corresponds to leaving out the term in parentheses, 
proportional to {C/CcurvY ®- However, this means that the energy losses of the 
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charged particle to curvature radiation are small compared to losses due to the term 
oc '-)d?'~^/dt^ in ([T]). These energy losses can also be associated with bremsstrahlung 
losses, due to the oscillatory motion of the particle. These will differ from zero only 
if the particle energy varies with time. In contrast, the curvature losses depend only 
on the radius of curvature of the particle trajectory and the particle energy. There is 
clearly no dependence of the intensity of curvature radiation on time derivatives of the 
energy of the charged particle. Thus, when Co ^ Ccurv, the particle energy decreases 
exponentially with a decay time constant Q, and the losses themselves occur mainly 
due to bremsstrahlung. 

Now let sq ^ 1. Obviously, in this case, the most interest is presented by the 
behavior of C on relatively short times t, since C decreases with time and will become 
smaller than Ccurv at sufficiently long times, so that the decay takes on an exponential 
character. At small times t, we have 



t 



e(t) = eo 1 (Co > t < 



where 



r, = 3. (11) 



^5 



We can see that, if Cq ^ Ccurv, there will be a linear decrease in the particle energy on 
times appreciably less than Tp. When t is comparable to r^, the dependence acquires the 
form 

/ t \ ^1/3 

e{t) = £0 (^1 + 3- J (Co > Ccurv. t^Tp). 

With further growth in t, the decay becomes a power law: 

^(^) = (3—1 (Co > Ccurv, Tp < t < Td). 

Here, the meaning of the time Tp becomes clear, as the time during which the decrease of 
the energy of the charged particle, which is initially linear, acquires a power-law character. 
On times of the order of r^, the dependence again changes: 

e{t) = [e'"^^ - ly'" (Co » Ccurv, t ~ r,). 

On still longer times, the dependence becomes exponential: 

e{t) = e-"^^ (Co > Ccurv, t > r,). (12) 

In contrast to the case (fTOl) considered earlier, when Co ^ Ccurv, the particle's energy 
losses are due primarily to curvature radiation. This process dominates on times t < Td 
and has a power-law character. When t > t^, the decrease in the power energy becomes 
exponential, indicating a transition to a regime in which the main contribution to the 
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energy losses is made by bremsstrahlung associated with the acceleration experienced by 
the particle during its oscillatory motion. It is clear that this transition itself occurs on 
time scales of the order of r^. However, we will refine this somewhat by using ([Tj) to find 
the time Tcurv for which C = Ccurv] i-e., when e{Tcurv) = 1: 

ln2 , . 

Tcurv = — Td ^ 0.23 Td. (13) 

Note that r^urv is completely independent of the initial particle energy Eq. This means 
that the charged particle, which initially has an energy ^ 1? loses a large fraction of this 
energy over the time Tcurv, even if its initial energy is very large. The particle will possess 
the energy e = 1 at time t = r^urv, after which there will occur the usual exponential 
decay with time constant r^. It is not surprising that the dependence e{t) [see f|T2|) ] also 
does not contain the initial particle energy at times t > r^. This is due to the fact that 
Tcurv is roughly a factor of four smaller than Td, and e becomes equal to unity over a time 
Tcurv, independent of Eq. Over long times, the relation e^Tcurv) = 1 begins to play the role 
of an initial condition for the further exponential decay, as is shown by ( [T2|) . 

All the results we have obtained by averaging ([T]) are valid when the adiabatic ap- 
proximation (|3]) is satisfied. As follows from ([6]), this condition yields an upper limit for 
the amplitude of the oscillations of the charged particle: 

A « A^ax = c.-^/ ^ 2.6 u-'/'p'/\ (14) 

The condition for the adiabatic approximation can also be interpreted as follows. The 
oscillation period should clearly be smaller than the characteristic decay time for the 
energy of the charged particle. At first glance, it seems that we should take the exponential 
decay time Td or the time scale T^urv for curvature losses as a characteristic time, but this 
is not so. We must take the more rapid time scale for the energy variations, which, in 
our case, is the time Tp for the transition to a power-law decay. Indeed, if we write the 
condition 

T < Tp, (15) 

the resulting upper limit on the oscillation amplitude coincides fully with the condition 
f|T^ . As a consequence, as a criterion for applicability of the adiabatic approximation, 
we can use any of conditions Q, (fT^ . or (|T5l) with equal success. 

Let us now estimate to order of magnitude all the quantities introduced above that 
characterize the time dependence of the energy of a charged particle oscillating near the 
force-free surface. We will first estimate the energy Ccurv'- 

Ccurv ~ {—f:, — ) • (16) 



For a characteristic surface magnetic field B ~ 0.01 — 0.1, neutron-star radius R ~ 10^^ 
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and light-cylinder radius Rl ~ 10^^—10^° (we will use these values for all our estimates 
below), Ccurv is, to order of magnitude, 10^. We can also introduce the corresponding 
amplitude 



which is Acurv ~ 10^^—10^^ (~ 1 m in dimensional units), for the same magnetic field, 
neutron-star radius, and light-cylinder radius. 

The time constant for the exponential decay. 



has characteristic values 10^^—10^^ (~ 10—10^ s in dimensional units). The time Tcurv 
has the same order of magnitude, since relation fll3p is satisfied. Note that the time 
exceeds the characteristic rotational periods P for pulsars. However, if the field becomes 
comparable to the critical fields, -B ~ 1, as is true for magnetars, the time can be even 
smaller than the rotational period of the neutron star. Note, by the way, that the ratio 
Trf/P depends only on the surface magnetic field 5, not on P. 

The order of magnitude of the time Tp cannot be fixed as well, since it depends on the 
initial particle energy Eq. However, we can easily find the range in which Tp lies. The upper 
limit of this range is given by the condition 1, oi equivalently Cq ^ Ccurv^ which is 
the only case it makes sense to introduce the concept of the time Tp. This means that, in 
any case, Tp <^ Td. We will determine a lower limit for the range of Tp values by taking 
the maximum possible value of the energy Cq. It is obvious that this value is clearly less 
than the maximum attainable Lorentz factor 70 that the particle will possess far from the 
force-free surface. Using the fact that 70 ~ 10*^ and assuming Eq ~ ^^j Ccurv ~ 10'', we find 
that the inequality Tp ^ lO'^^r^ is always satisfied. This lower limit is made with some 
reserve, since, in reality, the oscillating particle cannot possess an energy Ccurv ~ 10^: 
if the particle possessed such an energy far from the force-free surface, its energy after 
its capture would be somewhat lower, since the particle loses some of its energy in its 
approach to the force-free surface. 

Let us also find to order of magnitude an upper limit Amax for the oscillation ampli- 
tude A such that adiabatic approximation ([3]) is still applicable. Using (|T^ . we obtain 



Thus, Ainax ~ 10^^—10^^ (of the order of several tens of meters in dimensional units). 
This amplitude corresponds to the oscillation period Tmax ~ A, which corresponds to 
^max ~ 10^®— lO^'' s and v ~ 10—100 MHz in dimensional units. We can easily find the 




(17) 




(19) 
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maximum energy of the oscillations using the amplitude Ay, 

X 1/7 



Cmax ~ (^-^J , (20) 

which is Cmax ~ 3 X 10^ — 10''. This enables us to refine our lower limit for Tp. Taking as 
the maximum energy Eq the ratio Cmax/Ccurv ~ 10^ — 10'^, we obtain the more realistic 
lower limit Tp ^ lO^^r^. 

Finally, let us consider the energy losses of the charged particle in the case of non- 
relativistic oscillations, which is realized when the oscillation amplitude becomes less than 
Inro — l/uj- The non-relativistic Dirac-Lorentz equation has the form 

where the term containing the third time derivative of the coordinate is the well known 
expression for the force of radiative friction. The corresponding characteristic equation 
has three roots, one of which is real and positive and the other two of which are complex 
conjugates. The positive root must be thrown out, since it exhibits the effect of self- 
acceleration of the particle. The time dependence of the coordinate has the form 

x\ = Aoe"*/^""" cos cot, 

where Aq is the initial amplitude and the decay time constant is 



' nro o • 



As expected, to order of magnitude, this coincides with the decay time constant for 
the case of ultra-relativistic oscillations. 



3. CAPTURE OF CHARGED PARTICLES 

In jol, we found that a charged particle (electron or positron) far from the force-free 
surface moves in such a way that its Lorentz factor is determined only by its coordinates, 
and takes on values 70 ~ 10^. However, close to the force-free surface, the electric field 
becomes small and the time for variation of the particle's energy, tq = 3p^/8a7Q, becomes 
comparable to the characteristic distance (in dimensionless units) for variations in the 
field — in the case considered, the distance to the force-free surface. This distance Ic is 
determined in a self-consistent way from the relation tq = 1^, the longitudinal electric field 
E\\, on which tq depends, must be expressed in terms of the distance Ic, measured from 
the force-free surface along the magnetic-field line along which the particle moves. We 
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will use a linear approximation, assuming that Ic is sufficiently small that 

E\i{l) = -u^l, (21) 

where we have formally introduced = —dE\\/dl. Here, the derivative dE\\/dl must be 
taken at the point where the magnetic-field line along which the particle moves intersects 
the force-free surface. It is easy to see that uj is the frequency of non-relativistic oscillations 
of the particle. Further, we immediately obtain the capture length 

Note that the capture length coincides to within a multiplicative coefficient with the 
maximum oscillation amplitude Amax fll4p . determined by the adiabatic approximation 
condition. This suggests a transition from ultra-relativistic, quasi- stationary motion of 
the particle far from the force-free surface, determined by the energy-balance condition, 
to oscillatory motion of the particle near the force-free surface. Qualitative reasoning is 
sufficient to provide more insight into this. The capture length 1^ is nearly a factor of 
three smaller than the amplitude Amax- This suggests that a particle immediately makes 
a transition to the adiabatic oscillation regime (E]) upon its capture. To verify this, we will 
estimate the capture amplitude A^. — the maximum distance a particle deviates from the 
force-free surface when it passes through this surface. We suppose here that the charged 
particle was initially created far from the force-free surface, began moving toward this 
surface, and overshot the surface upon its first intersection with it. Let us consider a 
particle that has not yet intersected the force-free surface, located a distance Ic from this 
surface. We can estimate the Lorentz factor 7c of the particle at this point using (121 p and 
setting I = Ic'- 

At the same time, assuming the condition (|3]) is satisfied, we can calculate the first 
integral C, setting 7 = 7c and x'^ = Ic- 

We then have 

Ac = 3/c~2.9a;-^V^^- (22) 

The capture amplitude A^ essentially coincides with the maximum amplitude of adiabatic 
oscillations Amax (I14p . Consequently, we can take a charged particle that has passed 
through the force-free surface and traveled beyond it a distance Af. to be captured, after 
which it begins to oscillate in the adiabatic regime that we analyzed in detail above. 
For oscillations with an initial amplitude Ac, the time reaches its lower limit, ~ 
(10-9 - 10-^)r,. 
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For completeness, let us find the non-relativistic oscillation frequency u. Let us turn 
to ([2]), which determines w^. Calculation of the frequency requires computation of the 
scalar product B ■ V(E ■ B). In principle, there is no hindrance to calculating this product 
directly using the expression for the operator V in spherical coordinates. However, we 
will use a somewhat different approach that enables us to obtain the desired expression 
in a more compact form. We introduce the unit vectors = r/r, = m/m, and 
e„ = dem/ddmi where r is the radius vector, m the magnetic moment, and Om the angle 
between the magnetic axis and the rotational axis gq = O/fi of the neutron star. 
We can then determine the polar angles 9, 9', and 9" using the relations cos^ = ■ e^, 
cos 9' = er ■ e^n, and cos 9" = ■ e„ and the corresponding azimuthal angles and ip". 

Recall that the dipolar magnetic field is potential: B = — V$. Here, we have introduced 
the magnetic potential $ = mcos^'/r^, which is a harmonic function: A$ = 0. This 
enables us to write the desired scalar product as follows: 



B-V(E-B) = i 



$ A(E- B) - A(<l> E- B) 



The calculation itself is carried out by expanding the terms in the functions P™(cos 9) e*"*'''/r', 
where P^{cos9) are associated Legendre polynomials, and further calculating the Lapla- 
cian of each of these functions: 

A (^^Pr (cos 9) e^-^^ = ^Pr(cos 9) e^™^ [/(/ - 1) - n(n + 1)] . 

In place of the usual angles {9, (f), we can take the angles {9', (f') in a spherical coordinate 
system with the polar axis e.^, or the angles {9", ip") in a coordinate system with the 
polar axis e„. Without presenting the computations themselves, which are nevertheless 
cumbersome, we can write the final result as 



.3 



B- VfE-B) = ^cos^' 



^9 



, d2. d2 

COS 9" sin 9m + ^4 cos 9 (ll cos^ 9' + 3) 



(23) 



Substituting (!23l) into ([2]) yields for the square of the non-relativistic oscillation 
frequency 

2 km ( R\\ ^ ^,2cos2 0' + 3 , ^, 

w = — z— 4 COS 6' COS 6*- — -, (24) 



P3 \rffj 3cos2^' + l 

where r^fg is determined by the equation for the force-free surface 

cos 9 cos^ 9' 



r^,,=R' (1-4 



\ sin cos 

Here, we have used the relation for the square of the magnetic field = (3 cos^ 9' + 
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1) rn? /r^. Let us also present an expression for the radius of curvature of the field lines: 

_ r(l + 3cos^ e'f'^ 
^ ~ 3sin^'(l + cos2^')' 

If we must calculate p at some point on the force- free surface, we must set r = rjfs in 
this formula. 

As a rule, in the entire treatment above, we have not been concerned at all with what 
type of particle is oscillating near the force-free surface — an electron or a positron. We 
have assumed only that the sign of the charged particle is such that the force exerted on 
the particle by the longitudinal electric field Ey will act to return it toward the force-free 
surface. Only in this case is it possible to speak of oscillations of the particle. It is now 
straightforward to refine this question. All the expressions obtained above assumed that 
the oscillating particle had a positive charge. Consider formula flMj) for the non-relativistic 
oscillation frequency. When u"^ > 0, we have a positron oscillating with the frequency u. 
Accordingly, an electron cannot oscillate under these conditions, because the electric force 
will repel it from the force-free surface. If we formally suppose that u"^ < 0, then, on the 
contrary, oscillations of electrons will be occur. The non-relativistic oscillation frequency 
for an electron will then be a/^cJ^, and we can take the sign of the oscillating charged 
particle to coincide with the sign of o;^. 

As follows from f l24|) . the non-relativistic oscillation frequency vanishes if and only 
if 6 = 71 / 2 or 6' = 7t/2. The first equality corresponds to the equator of the neutron 
star, and the second to the star's magnetic equator. The equator and magnetic equator 
separate the force-free surface into regions in each of which oscillations of particles of 
only one sign can occur. The sign of the charged particles changes when passing through 
the line formed by the non-isolated equilibrium positions. Simultaneously, the equalities 
9 = 71/2 and 9' = 7i /2 define a line at each point of which a; = 0. However, the signs of the 
oscillating charges coincide on the open sheets of the force-free surface, which are divided 
by this line. If 6'^ < 7r/2, positrons will accumulate in the arched parts of the force-free 
surface, and electrons in its open sheets. Note that u"^ depends on the wavenumber fc, 
and therefore also on its sign. If > 7r/2, we must formally consider the case when 
the angle between the axes of an inclined rotator is n — 9m, also changing the sign of 
the angular velocity Q. Consequently, in this case, the signs of the charges will reverse, 
so that electrons accumulate in the arched parts and positrons in the open sheets of the 
force-free surface. 

4. PARTICLE TRAJECTORIES AT THE FORCE-FREE 

SURFACE 

The motion of a particle captured by the force-free surface is the sum of the drift motion 
of some driving center that moves exactly along the force-free surface and the oscillatory 
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motion of the particle about this driving center. Let us find the exact trajectories of the 
driving center in the rotating coordinate frame. 

Recall that the velocity of the driving center, Vq, can be represented as the sum of 
the drift velocity perpendicular to the magnetic field and the longitudinal velocity f 'jb. 
The velocity Vq itself always lies in a plane tangent to the force-free surface. However, 
this velocity possesses another important property. It is easy to see that 



Here, we have introduced the effective electric field E^'^-^ = E + v^^ x B, where E and 
B are the Deutsch electric and magnetic fields and wtr = f2 x r is the translational 
velocity. We can see that the velocity vector Vq of the driving center is always orthogonal 
to the effective electric field 'Ei'^^f . The reason for this is that, first, the drift velocity 
is orthogonal to the vector E"^-^-^ and, second, the driving center always remains on 
the force-free surface during its motion, Yff^ ■ B = 0, so that the longitudinal velocity 
component is also orthogonal E^-^-^. 

Let us consider the problem of finding an equation for the surface at each point of 
which the effective electric field Y/^^^ is normal to the surface. We will first consider the 
inverse problem. Suppose we have some two-dimensional surface specified by the equation 
^(r, 9, Lf) = const. Assuming that the function ^ is different iable, we introduce the vector 
field V^. Then, as is well known, this vector field is orthogonal to the surfaces of equal 
values of the function ^. In particular, the gradient at an arbitrary point lying on the 
specific surface ^ = const is orthogonal to this surface. Let us now return to the original 
problem. If there exists a function ^ such that 



the collection of surfaces that are orthogonal to the vector field E*^-^-^, are specified by sur- 
faces of equal values of the function ^. This function exists and is given by the expression 



As above, here, (r, 6, (p) are spherical coordinates, R the radius of the neutron star, k = 
1/Rl the wavenumber, m the modulus of the magnetic moment, 6^, the angle between 
the rotational and magnetic axes, and 6' the polar angle relative to the magnetic axis. 
The potential ^ is determined to within an arbitrary real constant, which we do not write 
here in the interest of brevity. Expression fj28|) itself was obtained via direct integration 
of (|27|) . Equation ( 127|) is most easily verified directly. In contrast to E^-^-^, in general, the 
proper electric field E is not a potential field, and also contains a vortex component. The 



v' ■ = 0. 



(26) 



^eff 



(27) 




(28) 
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equation of the surfaces orthogonal to the vector field E'^-^^-^ has the form 



e = ^, (29) 

where € = const is an arbitrary real number if (129|) has a solution. 

First let a particle be located at some point Xq of the force-free surface. This point 
can be substituted into the corresponding number = ^(xq), calculated using fl25]) . 
Since the velocity of the particle Vq lies in a plane tangent to the force-free surface and 
is orthogonal to the vector Eq-^-^ at the point Xq, the velocity lies in a plane tangent to 
the equipotential surface ,^ = £o- Thus, the velocity vector Vg is tangent to the curve 
intersecting the force-free surface and the equipotential surface ^ = (to at the point Xg. 
This means that the particle moves along this intersection, as can easily be verified by 
calculating the total time derivative of the potential ^. The potential itself is clearly 
time independent, d^/dt = 0, so that its total derivative has the form d^/dt = Vq ■ V^. 
Recalling that the effective electric field E^'^-^ is a potential field [see fl27|) ] and that the 
velocity is orthogonal to [see fl2^ ]. we immediately obtain d^/dt = 0. Consequently, 
the potential ^ is conserved along the trajectory of the particle — more precisely, along 
the trajectory of the driving center — and plays the role of an integral of the motion. 
We conclude that, as the driving center moves along the force- free surface, it also moves 
along the equipotential surface specified by the potential at the initial point Xq. As a 
consequence, the trajectory of the driving center on the force-free surface is determined 
by specifying the coordinate of only one point through which we construct the trajectory, 
and not necessarily the initial point. 

By virtue of the above, the task of finding the trajectory of the driving center reduces 
to finding the intersection of the force-free surface and the collection of equipotential sur- 
faces (12^ . We introduce the function ({0,ip) = C^{rffs{9,(p),9,ip)R/km, where rffs{9,ip) 
is the equation for the force- free surface and the potential ^ is determined by fl25]) . The 
desired trajectory is then given in implicit form by the equation 

ao,v) = €, (30) 



where, as above, (t is an arbitrary real number for which fl30|l has a solution. 

Figure 1 shows a phase portrait of trajectories in the coordinates {6, (f) on the surface 
in the rotating frame for the cases of an inclined and orthogonal rotator. Here and below, 
we assume that the azimuthal angle corresponding to the magnetic axis is zero. We 
can see that there are no extended trajectories having their beginning and end points 



on the magnetic equator. This differs from the results of Jackson [10|, who asserts that 
some trajectories on the force-free surface of an uncharged orthogonal rotator end at the 
magnetic equator, giving rise to a return drift motion at the stellar surface, which sweeps 
out charged particles. This could hinder the formation of a force-free magnetosphere. 
However, Jackson assumed that the velocity of a particle at some point of the force-free 
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surface is equal to the projection of the drift velocity onto the tangent plane at that point. 
We have elucidated above that this is not the case. 

Figure 2 presents a three-dimensional image of the trajectories of the driving center 
on the force-free surface for the case of an inchned rotator. We introduce for convenience 
in our further analysis the potentials of the isolated equilibrium points, 

l (lTcosg,„)-V^ 
^3(3±cos^JV2 

and the surface potential 

2 

Cr = g cos 61^. 

At the cupola-like parts of the force-free surface, the trajectories have the form of closed 
cycles surrounding the equilibrium position. The equilibrium position itself corresponds 
to the potential C_, and the collection of cycles to the interval of the potential < C < C-- 
The equilibrium positions at the equator and magnetic equator are characterized by the 
surface potential (r. Two regions can be distinguished in the phase portrait at each of 
the two open sheets. The first region is analogous to the region of the arches, and also 
consists of closed cycles surrounding the equilibrium position. The equilibrium position 
corresponds to the potential C+ and the collection of cycles to the potential interval 
C_i_ < < 0. However, the cycles themselves become more elongated and tend to go 
out to infinity as ( approaches zero. The scparatrix corresponds to the potential C = 0- 
The scparatrix itself represents an open trajectory that goes out to infinity, but also 
encompasses the equilibrium position. The second region in each of the open sheets is 
formed of a collection of trajectories corresponding to the potentials < C < (^r. None of 
these trajectories lies entirely within a single sheet. As Fig. 2 shows, all the trajectories 
pass from one sheet to the other, intersecting the line 9 = 9' = tt /2 and closing around the 
neutron star. As ( approaches (r, these trajectories come closer to the magnetic equator. 

The direction of the trajectories is not difficult to determine. For this, we present 
the components of the velocity v^: 



'■'±r — ^ ^rb(p 



2\ r^J' 



'±9 



Qr 



r 



cos 9br + (1 ^ ) sin 9 be 



(l - — ) sin 9 (^-bl + bi^ + cos 6 bM 



where br, be, and b^ are the components of the unit vector b. By virtue of the structure of 
the trajectories, we can consider values of the velocity only at the points of intersection 
of the force-free surface and the plane cp — {0, tt}. Here, b^ — 0, so that v'j_j. — v'j_0 — 0. 
At the points considered, the velocity of the driving center fully coincides with v^, so 
that I'll = 0. To determine the required sign of the component v'j_^, we turn to the 
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trajectories that are infinitely close to the magnetic equator, cos^' = 0. The expression 
1 — i?^/r^ is quadratic in the small quantity cos 6*', while br is linear in this quantity. 
Consequently, the component v'j_^ changes its sign at the intersection with the magnetic 
equator, with sgn{v'j_^) = — sgn{br cos 6) in a small vicinity of the magnetic equator. 
Knowledge of this fact is all that is required to reconstruct the direction of the trajectories 
everywhere on the force-free surface using the principle of continuity. We have at the 
points of intersection of the force-free surface and the = half-plane v'j_^ > when 
9 e (0, 9^m/2) U (7r/2 + 6^/2, 7i/2 + Om) and v'^^ < when 9 E (6^/2, ft^) U {n/2, 7r/2 + 
9m/2) U (7r/2 + 9m, vr). We have at the points of intersection of the force-free surface and 
the ip = 71 half-plane v'j_^ > when 9 G {11/2 — 9^,7^/2 — 9^/2) U (vr — 6*^/2, vr) and v'j_^ < 
when 9 G (0, 7r/2 - 9rn) U (7r/2 - 9m/2, tt /2) U (tt - 9^, tt - 9m/2). 

Thus, in the half-space cos 6' > in the arched part of the force-free surface and on 
the open sheet inside the C = separatrix, the motion along the trajectories is clockwise 
as seen by an external observer (of course, the line of sight passing through the force- 
free surface is directed along the radius toward the center of the star). In the half-space 
cos 9 < in the arched part of the force-free surface and on the open sheet inside the 
separatrix, the motion along the trajectories is counter-clockwise as seen by an external 
observer. On the open sheets outside the separatrix (for trajectories intersecting the line 
9 = 9' = 7r/2), the motion passes from the half-space sinyj < to the half-space sinyj > 
when cos^ > 0, and from the half-space simp > to the half-space simp < when 
cos 9 < 0. All the above refers to the case 9m < 7r/2. If the inclination angle of the 
rotator is x > 7r/2, we must treat the rotator as having an inclination 9^ = tc — x, with 
its polar axis, from which all polar angles are measured, being directed opposite to f2, 
and change the sign of Q. In Fig. 2, this corresponds to changing the direction of fl. It 
is obvious that the directions of all trajectories will be reversed. 

Thus, all trajectories of the driving center on the force-free surface are closed and 
lie in a finite region, apart from the separatrix C = 0. As a charged particle moves along 
its trajectory, the electromagnetic field, and therefore the parameters of the particle's 
oscillation, change. However, the driving center moves with a drift velocity ~ 10~^, and, 
as follows from f|T9|) . the oscillation frequency exceeds 10 MHz, even for ultra-relativistic 
oscillations, so that the particle traverses a distance of the order of several millimeters, 
appreciably less than R, over the oscillation period. As a consequence, the variation of 
the parameters of the particle oscillation during its motion along the force-free surface is 
adiabatic. Note that the sign of u"^ does not change in the motion along the trajectory 
[see the discussion following (1251) ]. This means that, if a particle falls into the force- free 
surface, it is not able to leave, since, first, the trajectory of its driving center is closed and, 
second, instability of the particle oscillations along the magnetic field does not develop 
due to the constancy of the sign of cu^. 
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5. ACCUMULATION OF ELECTRON-POSITRON 

PLASMA 



Let us estimate the rate of accumulation of electron-positron plasma at the force-free sur- 
face. We will consider the situation when the charged particles — electrons and positrons — 
move in a specified electromagnetic field, playing the role of test charges. This is possible 
only while the density of charges in the magnetosphere is appreciably less than the Gold- 
reich- Julian density pcj- However, the rate at which pairs are created in the magneto- 
sphere is constant, so that, with time, particles are accumulated at the force-free surface. 
Consequently, a region appears in the vacuum magnetosphere near the force-free surface 
where the charge density becomes comparable to pcj over some finite time. As a result, 
the entire magnetosphere becomes filled with plasma. 

We first consider the initial stage in the restructuring of the magnetosphere from a 
vacuum state to a state in which it is filled with plasma. We choose an arbitrary point fq 
on the force- free surface and consider a small area on this surface containing this point. 
The size of this area is assumed to be small compared to i?, so that the area is virtually 
flat. With time, charged particles accumulate in this area, forming a symmetrical charge 
layer. 

We will consider the accumulation of plasma only for fairly short times, when the 
thickness of the formed charge layer is appreciably smaller than R. In this case, the plasma 
density differs appreciably from zero only near the force-free surface, so that there is not 
yet a global restructuring of the magnetosphere. Consequently, the outer electromagnetic 
field where the charged layer is located does not change. Let us find the proper electric 
field of the layer at some distance z from the force-free surface: 



where Pe{z') is the volume charge density a distance z' from the force-free surface and the 
distance z does not exceed H — the distance from the force-free surface to the boundary 
layer. Here, it is important that the charge layer is symmetrical, pe{—z') = Pe{z')- Pro- 
viding an equilibrium state for such a layer requires that each charged particle belonging 
to the layer not be accelerated along the magnetic-field lines. In other words, the outer 
longitudinal electric field fl2T|) must be fully balanced by the projection of the electric field 
of the layer (ETj) onto the direction of the magnetic field. Introducing the angle if) between 
the normal to the force-free surface and the direction of the magnetic field. 



z 




(31) 







B- V(E-B) 
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The relation 

I cos l/) 



47racosV' J Pe{z') dz = uP'l (32) 







is then satisfied for arbitrary /. Differentiating both sides of (132|) with respect to /, we 
obtain 



= I TT- (33) 

Ana cos^ ip 

Thus, the charge density pe in the layer does not depend on the distance from the 
force-free surface. Recall that we are measuring the charge density in units of e/A'^. This 
is determined only by the square of the non-relativistic oscillation frequency and the 
angle ip, whose values are taken at the point tq lying in the small area of the force-free 
surface we are considering. The sign of the charge density pe coincides with the sign of 
w^, in agreement with our conclusions following (125 p . In general, pe does not coincide 
with pgj- This flows from the fact that, if the charge density in the layer were everywhere 
equal to the Goldreich- Julian charge density, all the particles in the layer would be fully 
co-rotating, so that they would be stationary in the rotating frame. In the layer, only the 
longitudinal electric fleld is zero, while the flelds perpendicular component differs from 
zero, because the electric fleld of the layer does not necessarily cancel out the perpendicular 
component of the external fleld. For example, because the electric fleld of the charged 
layer is equal to zero at the point tq, the total electric fleld has only a perpendicular 
component, which gives rise to a drift of the charged particles along the force- free surface, 
along the trajectories depicted in Fig. 2. This means that the forming charge layer is 
made up of differentially flowing currents on the force-free surface. 

Thus, a charge layer forms near the force-free surface with time, whose thickness 
if is a function of the point tq and the time t that has passed since the start of the 
flUing of the magnetosphere with plasma. It is important that only charges of the same 
sign as the charge density pe contribute to the value of p^ at each point of the force-free 
surface. Because of this, the layer is fully charge-separated in the initial stage of flUing 
of the magnetosphere, in accordance with our conclusions following (125|) . Consequently, 
the number density of the particles ne coincides with the plasma density: Ue = \pe\- 
At each point of the force-free surface, we introduce the current density of the particles 
j = dN/dSdt, which coincides with the electric-current density in dimensionless units. 
The equation used to flnd the thickness of the charge layer then has the form 

Equation can be used to estimate the time to flU the entire vacuum magneto- 
sphere with electron-positron plasma. When the thickness of the charge layer H becomes 
comparable to the characteristic size of the inner magnetosphere, i.e., to the radius of the 
star R, there is a substantial restructuring of the entire magnetosphere. The flUing time 
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is then 

Rncj 

= — — • 

J 

Here, we have taken into account that the plasma density in the layer is of the order of 
the Goldreich- Julian density (ng ~ ncj), as follows from fl33p . The current density of the 
particles incident onto the force-free surface j is determined by the mechanism creating 
pairs in the magnetosphere. Gamma-rays capable of producing electron-positron pairs in 
the magnetosphere magnetic field can both be incident from outside, as cosmic rays, and 
be produced from soft thermal photons radiated by the stellar surface that are Compton 
scattered by energetic particles. The main question here is how many pairs are produced 
by a single photon with energy exceeding 1 MeV. Since the created electrons and positrons 
are very rapidly accelerated to energies 70, there is a chain multiplication of pairs. The 
number of created pairs per primary photon depends exponentially on the characteristic 
size of the inner magnetosphere: exp /x, where fi = R/It.. Here, /s is the total mean free 
path of an energetic particle, including both curvature radiation and pair creation. 

Studies of processes associated with the stationary creation of plasma in the polar 



magnetosphere Uj indicate that the mean-free path Zs is of the order of 100 m. This 
means that the multiplication coefficient can reach gigantic values of the order of exp 100 ~ 
10^^. However, in view of the presence of the exponential factor, resolving this question for 
the case of a vacuum magnetosphere requires a more careful analysis, which falls outside 
the framework of this paper. 



6. CONCLUSION 

We have shown that as a charged particle that is initially accelerated to a Lorentz factor 
7o ~ 10^ approaches the force-free surface, the quasi-stationary condition is violated, and 
the particle passes through the force-free surface with a Lorentz factor Cmax ~ 10^ — 10^ 
fl2U]) . After passing through the force-free surface, the particle deviates from this surface 
by a distance Ac ~ Amax ~ 10 — 100 m (fT^ . (fT9|) . (122|) . Further, the particle undergoes 
adiabatic ultra-relativistic oscillations at a frequency z/ ~ 10 — 100 MHz ( 1T9|) . These 
oscillations decay due to energy losses to radiation, and their frequency grows. The 
energy of the oscillations falls first linearly, then, over a time Tp ~ 10~^ — 10"^ s [see fllip 
and the comments following f l20|) and fl22l) ]. according to a power-law decay, due to losses 
to curvature radiation. Over times of the order of Tcurv ~ ~ 10 — 1000 s ffT^ . f lTS]) . 
the decay becomes exponential, with the time constant , with the main contribution 
to the energy losses now made by bremsstrahlung. During the change in the decay regime, 
the particle possesses a Lorentz factor Ccurv ~ 10"^ ([8]), (|T6l) . and the oscillation amplitude 
is Acurv ~ 1 m fll7p . Further, the ultra-relativistic oscillations of the particle continue to 
exponentially decay. When the amplitude Inro ~ 1 cm is reached, the oscillations become 
non-relativistic and harmonic, and their frequency is z/ ~ 1 — 10 GHz. Simultaneous with 
this oscillatory motion, the particle undergoes regular motion along the force-free surface 
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(Fig. 2). As a rule, the trajectory of the driving center is closed in the rotating frame, 
and its velocity is of the order of the drift velocity. 

The gradual accumulation of particles at the force-free surface leads to the formation 
of a fully charge-separated layer of electron-positron plasma. The volume density of this 
plasma depends on the coordinates (6*, ip) on the force-free surface, but not on the layer 
thickness. This layer has a sharp boundary, and the plasma density in the layer is com- 
parable to the Goldreich- Julian density. In the case of a constant rate of creation of pairs 
in the magnetosphere, the layer thickness grows monotonically. Over a finite time, the 
layer thickness becomes comparable to the radius of the neutron star, accompanied by a 
restructuring of the magnetosphere. This is the characteristic time for filling of the mag- 
netosphere with plasma. The efficient acceleration of the particles gives rise to a chain of 
creations of secondary pairs. The resulting high multiplication coefficient for pair creation 
may be compensated for by the smallness of the flux of cosmic-ray photons. According to 
the observations [12I, [3], the diffuse background of Galactic photons jph having energies 
exceeding 1 MeV is of the order of 10~^ cm~^s^^. Thus, the ratio jph/cncj — 10^^^ is 
very small. However, this does not hinder the rapid filling of a pulsar magnetosphere with 
electron-positron plasma if the pair-creation multiplication coefficient is sufficiently high. 
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Fig. 1. Phase portrait of the trajectories of charged particles at the force-free surface in 
the coordinates {6, if) in the rotating frame for rotators whose rotational and magnetic 
axes are mutually inclined by 7r/3 (left) and 7r/2 (an orthogonal rotation, right). 
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Fig. 2. Trajectories of charged particles on tiie force-free surface in the rotating frame for 
the case of a rotator with its rotational and magnetic axes inclined by 7r/3. 
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